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1. Introduction 



While the work of Gu and Zhu [10] establishes the existence of Ricci flow solutions 
that form Type-II singularities, it tells us very little regarding the details of the 
evolving geometries of such solutions. The numerical work of Garfinkle with one 
of the authors [8, 9], together with the formal matched asymptotics derived by 
all three authors [2], strongly suggest what some of these details might be, at 
least for those solutions (§ n+1 , g(t)) that are rotationally symmetric and involve a 
degenerate neckpinch. However, these works do not prove that solutions with the 
prescribed behavior exist. In the present work, we prove that indeed, for each of the 
prescribed models of the evolving geometry near the singularity discussed in our 
work on matched asymptotics [2] , there is a Ricci flow solution that asymptotically 
approaches this model. Specifically, we prove the main conjecture stated in our 
earlier work [2, §7]. It follows from the results proved here and stated below that 
in every dimension n + 1 > 3 and for each integer k > 3, there is a Ricci flow 
solution that develops a degenerate neckpinch singularity with the characteristic 
profile (locally) of a Bryant soliton, and with the rate of curvature blowup given by 



for t approaching the singularity time T < oo, and for some constant C . 

In determining that there exist Ricci flow solutions that form Type-II singu- 
larities, Gu and Zhu [10] show that there are solutions in which the curvature of 
the evolving metrics satisfies limsup t y. T {(T — i) sup^gg^+i |Rm(x, = oo. Their 
work does not, however, determine the specific rate of curvature blowup in these 
solutions. Enders [7] has defined Type- A singularities as those with curvature blow- 
ing up at the rate (T — i)~ r , with r £ [1, |), and has questioned whether there are 
any compact solutions outside this class. Our work shows that indeed there are: as 
noted above, we prove that there exist solutions with curvature blowing up at the 
rate (T — t)~ 2+2 / k for all positive integers k > 3; all but the k = 3 solutions fall 
outside Enders' Type-A class of solutions. 

Noting the existence of solutions with these discrete curvature blowup rates, we 
are led to ask if there exist (compact) solutions — degenerate neckpinch or otherwise 
- that exhibit blowup rates other than these discrete values. Our work does not 
shed light on this question. We do note that for two-dimensional Ricci flow, in 
which case the flow is conformal and the conformal factor evolves by logarithmic 
fast diffusion, u t = A log it, Daskalopoulos and Hamilton [6] prove that there exist 
complete noncompact solutions that form singularities at the rate (T — t)~ 2 . 

The class of metrics we consider in this work are SO(n + l)-invariant Ricmann- 
ian metrics on the sphere § n+1 . We work with Ricci flow solutions g(t) for such 
initial data, and focus on solutions that develop a singularity at one or both of 
the poles at finite time T (proving that such solutions do exist). To define what it 
means for a singularity to be a neckpinch, we recall that a sequence {(xj, tj)}°^ 
of points and times in a Ricci flow solution is called a blow-up sequence at time 
T if tj /* T and if \ Rm(x j,tj)\ — > oo as j — > oo; such a sequence has a corre- 
sponding pointed singularity model if the sequence of parabolic dilation metrics 
gj(x, t) := | Rm(xj, tj)\ g(x, tj + | Rm(ij, tj)] -1 1) has a complete smooth limit. We 
say that a Ricci flow solution develops a neckpinch singularity at time T if there is 
some blow-up sequence at T whose corresponding pointed singularity model exists 
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and is given by the self-similar ly shrinking Ricci soliton on the cylinder 1 x §". 
We call a neckpinch singularity nondegenerate if every pointed singularity model 
of any blowup sequence corresponding to T is a cylindrical solution, and we call it 
degenerate if there is at least one blowup sequence at T with a pointed singularity 
model that is not a cylindrical solution. 

Rotationally symmetric nondegenerate neckpinches have been studied exten- 
sively by Simon [12] and by two of the authors [3, 4]. In the latter works, it is 
shown that there is an open set of (rotationally symmetric) compact initial mani- 
folds whose Ricci flows develop nondegenerate neckpinch singularities, all of which 
are Type-I in the sense that limsup (/ * T {(T — t) sup^ggn+i | Rm(i,f)|} < oo. Fur- 
ther, in the presence of reflection symmetry, [4] provides a detailed set of models 
for the asymptotic behavior of the geometry near a developing nondegenerate neck- 
pinch, with those models collectively serving as attractors for these flows. 

Unlike the formation of nondegenerate neckpinches, the formation of degenerate 
neckpinches in Ricci flow is expected to be an unstable property. This is evident in 
[8, 9] as well as in [10]: in all three of these works, one studies flows that develop 
degenerate neckpinches by considering one-parameter families of initial data such 
that for all values of the parameter above a threshold value, the Ricci flow solutions 
develop nondegenerate neckpinches, while for all parameter values below that value, 
there is no neckpinch singularity. The flows with initial data at the threshold value 
of the parameter are the ones that develop degenerate neckpinches. This instability 
leads us to use the somewhat indirect Wazewski retraction method [13] to explore 
the asymptotic behavior of the geometry near these degenerate neckpinches as they 
form. We discuss this in detail below. 

The matched asymptotics derived in [2] rely heavily on the imposition of a series 
of Ansatz conditions to characterize the formal solutions of interest. (It is through 
these Ansatz conditions that one builds the formation of degenerate neckpinches 
into the formal solutions.) By contrast, no such a priori assumptions are needed 
(or used) in the present work. Rather, having determined in [2] the nature and 
the explicit approximate forms (in regions near the degenerate neckpinch) for our 
formal solutions, we show here without imposing any further assumptions that there 
exist Ricci flow solutions that approach each of the formal solutions. It follows that 
degenerate neckpinches form in these Ricci flow solutions. 

Each of the solutions we consider here is characterized by an integer k > 3. 
As our results show, there is at least a one-parameter family of solutions corre- 
sponding to each value of k. In fact, our construction in Section 7 below reveals 
that degenerate neckpinches form in solutions starting from a set of initial data 
of codimension-fc in the space of SO(n + l)-invariant solutions. Besides determin- 
ing the rate of curvature blowup, the integer k also characterizes to an extent the 
detailed asymptotic behavior of the solution in a neighborhood of the singularity. 
While we discuss the details of this characterization below (in Section 3), we note 
here one important feature that depends only on the parity of k. If k is even, then 
the solution is reflection symmetric across the equator, neckpinch singularities oc- 
cur simultaneously at both poles, and the volume of (§™ +1 , g(t)) approaches zero 
at the time of the singularity. If on the other hand k is odd, then the neckpinch 
occurs at one pole only, and the volume of (S™ +1 , g(t)) remains positive at the time 
of the singularity. We note that in either case, for t < T, the curvature has local 
maxima both at one of the poles and at a nearby latitude sphere S™ where the neck 
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is maximally pinched. As t approaches the singular time T, the distance between 
the neckpinch sphere and the pole approaches zero, and the curvature becomes 
infinite simultaneously (albeit at different rates) both at the neckpinch sphere and 
at the pole. 

A detailed statement of our main results depends on the details of the asymptotic 
behavior of the formal model solutions. Referring to our discussion of this behavior 
below in Section 3, we can state our main theorem as follows: 

Main Theorem. For every integer k > 3 and every real number < 0, there exist 
rotationally symmetric Ricci flow solutions (S™ +1 , g(t)) in each dimension n + l > 3 
that develop degenerate neckpinch singularities at T < oo. For each choice ofn,k, 
and bk, the corresponding solutions have distinct asymptotic behavior. 
In each case, the singularity is Type-II — slowly forming — with 

jup i |Rm(x,t)|= (T _ t)2 _ 2/& 

attained at a pole, where C — C(n, k,bk)- 

Rescaling a solution corresponding to {n, k, bk} so that the distance from the 
pole dilates at the Type-II rate (T — i)^^ 1 ^ 1 /^, one finds that the metric converges 
uniformly on intervals of order (T — to the steady Bryant soliton. 

Rescaling any solution so that the distance from the smallest neck dilates at the 
parabolic rate (T— t)^ 1 / 2 , one finds that the metric converges uniformly on intervals 
of order y/T — t to the shrinking cylinder soliton. 

Furthermore, the solutions exhibit the precise asymptotic behavior summarized 
in Section 3, and they satisfy the estimates summarized in Section 4-4- 

Since the formal model solutions play a major role in this work, after setting up 
the needed coordinates and metric representations for our analysis in Section 2, we 
carefully review the nature of these formal solutions and their matched asymptotic 
expansions in considerable detail in Section 3. While the formal solutions we discuss 
here are the same as those analyzed in [2], we note that here we use somewhat 
altered coordinate representations in certain regions near the pole. In Section 4, 
we describe the general structure of the proof of the main theorem and outline its 
key steps. The technical work to carry out these steps is detailed in Sections 5-7. 

2. Coordinates for the four regions 

As noted above, in this paper we study SO(n + l)-invariant metrics g evolving 
by Ricci flow on § n+1 x [0, To), where To = To (go) € (0, oo] for initial data go. Each 
such metric may be identified with functions tp, ip : (—1, 1) x [0, To) — > R+ via 

(2.1) g(x, t) = tp 2 (x, t) (dx) 2 + ip 2 {x, t) 5can , 

where g can is the canonical round unit-radius metric on Smoothness at the 
poles requires that tp, ip satisfy the boundary condition (2.4) given below. Under 
Ricci flow, the quantities tp and ip evolve by 

- f ^ xx - Vx ^ x \ 
< ^ t V tpip p 2 ip ) ' 

4>xx <Px4>x , , , s tpl n-l 



DEGENERATE NECKPINCHES IN RICCI FLOW 



5 



respectively. This system is only weakly parabolic, reflecting its invariance under 
the full diffeomorphism group. Below, we remedy this by suitable choices of gauge. 

We are interested in proving the existence of solutions that are close to the 
formal solutions constructed in [2]; accordingly, we follow the notation of [2] in 
large measure. That paper describes solutions in four regions, the outer, parabolic, 
intermediate, and tip, using either of two coordinate systems. We do the same here, 
although we choose the alternate coordinate system to describe the intermediate 
and outer regions. For clarity, we review both systems below. 



2.1. Coordinates for the parabolic region. We call the SO(n+l)-orbit {0} xS r ' 
the "equator" and denote the signed metric distance from it by 



s(x,t) := / ip(x,t)dx. 
Jo 

Then the metric (2.1) may be written as 
(2.2) g = {ds) 2 + ^ 2 {s{x,t),t)g c[ 

where d denotes the spatial differential. We write 

d 



d_ 
dt 



and 



dt 



to indicate time derivatives taken with x and s held fixed, respectively. With this 
convention, one has the commutator 

4> ss d 



' d 


d ' 


dt 





and the relation 



dt 



dt 



+ nl[ip] 



ip ds 
d 



ds : 



where I[ip] is the nonlocal term 

I[1>](8,t) 



ds. 



In terms of these coordinates, the evolution of the metric (2.2) by Ricci flow is 
determined by the scalar equation 



(2.3) 



dtp 
~di 



ipss - (n - 1) 



1> 



Smoothness at the poles requires that ip satisfy the boundary conditions 
(2-4) lk| x =±i=Tl. 

The quantity ip, which is effectively suppressed in these coordinates, evolves by 



d(logy) 
dt 



tpss 
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2.2. Coordinates for the other regions. We employ a different coordinate sys- 
tem in the other regions. In any region where tp s ^ (e.g., in the neighborhood of 
the north pole x = 1 where if> s < 0) one may use ip as a coordinate, writing 

(2.5) g = z(ip,t)- 1 (dij;) 2 + ip 2 g can , 

where 

z(ij{s,t),t) :=^ s {s,t)- 
The evolution of the metric (2.5) is determined by 

dz 



^ at 



where E^p is the quasilinear (but purely local) operator 

(2.7) := zzw -\zl + {n-l-z) Z ^ + 2(n - 1)^=^- 

We can split the operator £ into a linear and a quadratic component, 

£^[z] = £^[z] + Q^[z], 

where 

(2.8a) ^ M:=(n _i){^ +2 _L} ; 

(2.8b) Q^z] := zz^ _ _ ^ _ 2 (n - 1) * 

respectively. The quadratic part also defines a symmetric bilinear operator, 

(2.9) Q^[2i,2 2 ] := ^{^lO^vv + ^2(^1)^ - («i)v(^2)^} 

1 ,- 1 . -^2 

2 V ^ 2 ' 
In terms of this notation, one has Q^[z\ = Q^[z, z]. 

3. The formal solution revisited 

In [2], we present a complete formal matched asymptotic treatment of a class 
of rotationally symmetric Ricci flow formal solutions that form degenerate neck- 
pinches. These formal solutions serve as the approximate models which the solu- 
tions we discuss here asymptotically approach. Since we find it useful in our present 
analysis to work with different coordinate representations (namely z[u, r)) in the 
intermediate region than we use for the same region in [2] (effectively u(a,t)), we 
now briefly review some of the analysis of [2]. 

3.1. Approximate solutions in the parabolic region. Roughly speaking, the 
parabolic region is that portion of the manifold, away from the tip, where the 
geometry approaches a shrinking cylinder, and where the diameter of the neck 
has at least one local minimum. (We give a precise definition of this region in 
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equation (4.6) in Section 4.4 below.) As discussed in [2], it is useful in this region 
to work with coordinates consistent with a parabolic cylindrical blowup: 1 

(3 - 1) ° := v^-w-ty ° := 7hi- '■—**?-(>■ 

In terms of these coordinates, the Ricci flow evolution equation becomes 2 
with 

(3.3) iMfar) = U ^[' J ^ ) da. 



{^ + nI[u])u <J + \(u- l -) + {n 



o u{a,r) 

Setting u = 1 + v and linearizing at u = 1, we are led to 

dv 



— v a + v + { nonlinear terms }. 



This form suggests writing the solution to this equation by expanding v in Hermite 
polynomials h m , which are eigenfunctions of the linear operator 

(3-4) Ai= *-!.° +1 

do 1 2 do 

satisfying (A + fx m )h m = 0, where 

(3.5) Hm ■= - 1- 
Writing v in this fashion leads to the approximation 

oo 

(3.6) uw 1+ ^b m e-^ T h m {a). 

m=0 

This expansion can at best be an approximation to the actual solution, if only 
because the variable a is bounded. Note that one characterization of the parabolic 
region is that it is the (time-dependent) range of a for which the series terms in 
equation (3.6) are sufficiently small. 

We are interested in solutions for which the term with m = k is dominant for 
some specified k > 3. (Compare Ansatz Condition 2 of [2].) Thus the approxima- 
tion above takes the form 

u = l + b k e-^ T h k (a) + -- - , 

where the sign b k < reflects the fact that the singularity forms at the north pole. 

We normalize so that the leading term in hk{o~) is a k . We show below that, as 
they evolve, solutions of interest become C 1 -close to the formal solution 

(3.7) u = 1 + b k e-^ T a k + ■■■ 

in that portion of the parabolic region with \o~\ 3> 1. Assuming this for now, 
we compute an approximation of the formal solution in terms of the quantity z 
introduced above. To approximate z, we differentiate the expression for u with 
respect to cr, leading to an expression for z in terms of a. Using the relation 



4n Section 2, To > denotes the maximal existence time of a solution with initial data go- 
Here, T > denotes the singularity time of the formal solution we construct below, following [2]. 
In what follows, T is fixed, albeit arbitrary and unspecified. 

2 See equation (3.4) of [2]. 
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between u and a, we then write z in terms of u. The details of the calculation are 
as follows. 

Since tp = u^j2(n — 1)(T — t) and since s = o\jT — i, we have 
z = = 2(n - l)ul « 2(n - l)*^"^) V s " 2 
so long as \a\ ^> 1. On the other hand, we also have 

where 

(3 - 8) 7fc:= T = 2-fc- 

Thus for \a\ ^> 1, we get 

z^2(n~\)k 2 b\e-^-^a 2k - 2 

«2(»-l)^e-^ ({^f} lA e^) : 

= 2(n - l)fc 2 (-fe fe ) 2 / fe e- 2 ^ r (l - W ) 2 - 2 / fc . 

Hence at the interface between the parabolic and intermediate regions, where u is 
slightly smaller than 1, one has 

(3.9) 2«c fc e- 2 ^ r (l-u) 2 - 2 / fc , 
with Cfe defined by 

(3.10) c fe :=2(n~l)k 2 (-b k )V k . 

3.2. Approximate solutions in the intermediate region. The intermediate 
region is a time-dependent subset of the neighborhood of the north pole where 
— 1 < u a < and < u < 1. (We provide a precise definition in (5.12) below.) By 
equation (2.6), we know that 

dz 

Since V = \/2(ra - l)(T-t)u = i/2(n - l)e" r/2 u, we obtain 



(3.11) | 



1 P M_i 



2(ti- 1) J 2 

Noting the expansion (3.9) for z near u = 1, our first impulse is to look for 
approximate solutions of the form z rs e~ 27fcT Zi(u). Wishing to refine this approx- 
imation, with the goal of constructing lower and upper barriers Z- < z < z+ for 
the intermediate region, we are led by past experience to construct a formal Taylor 
expansion in time, namely 

(3.12) z = e-* 1 * T Zi(«) + e-^ kT Z 2 (u) + ■■ ■ = ^ e~ 2m ~< kT Z m (u). 

m>l 

We substitute this expansion into (3.11) and split E u [z] into linear and quadratic 
parts, as in (2.8). By comparing the coefficients of e" 2m7fcT in the resulting equation, 
we find that Z m must satisfy 

1 dZ 

-2mj k Z m = n , ' ,x { &u[Zm] + } ' Qu[Zi, Z m -i] } - -u- 
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This leads to the family of ODE 

(3.13) \( u l ~ U ^~^T + I + 2m 7k) Z m = ~ 2Cn- 11 ^ Zm-i], 

from which the Z m can be computed recursively, up to a solution of the associated 
homogeneous equation (obtained by replacing the RHS of equation (3.13) with zero). 
The general solution of the m th homogeneous equation is 

(3.14) Z m . hom (u) = c m u- 2 (l - u 2 ) 1+2mi \ 

where c m is arbitrary. 

For m = 1, the RHS of equation (3.13) vanishes. Therefore Z\ satisfies the 
homogeneous equation, whereupon setting c\ — c k yields 

(3.15) Z x (u) = Z 1Mom {u) = c k u- 2 (l - u 2 ) 1+2 ^ k . 

The function Z 2 is harder to compute (we have not tried to find a solution in closed 
form). However, the expression (3.15) for Z\ shows that Z\ « Cfc2 1+27fc (1 — u) 1+2 ~ lk , 
Zi « -c fe (l + 2 7fe )2 1+2 ^(l - u) 2 ^, and Z/' » c k { lk + 2 7 2 )2 2+27 *(l - w) 27 *" 1 as 
u /• 1. Thus we calculate 



QuIZlZx] = Q u [^] « -c 2 2 1+4 ^(l - 4 7 2 )(1 - u ) 4 ^. 

It follows that for some constant c, the solution Z 2 of equation (3.13) takes the 
form 

(3.16) Z 2 (u) = (c + o(l))(l - u) 47fc as u /* 1. 

The solution .Z^.hom of the homogeneous equation is dominated by (1 — u) 1+47fc , 
hence is smaller than the RHS of (3.16) for u /• 1. Therefore, all solutions i?2 of 
equation (3.14) satisfy (3.16). 

It follows from these calculations that, near the neck, where u ~ 1, the first two 
terms in our formal expansion of z predict that 

(3.17) z = c k e~ 2 ^ r u- 2 (\ - u 2 ) 1+27fc + 0(e- 47fcT (l - uf lk ). 

Since "f k < |, we have 1 + 2 7 j. > 4 7 j.. Thus for fixed r, the remainder term 
becomes larger than the leading term as u / 1. In [2], we characterize the in- 
termediate region using the coordinate system (2.2); here, we identify it in terms 
of the system (2.5) as the region where our formal solution e~ 2lkT Z\ is a good 
approximation to the true solution. The formal solution is a good approxima- 
tion only so long as the remainder is smaller than the leading term. This is 
true where e _47fer (l — u) 4 "' k <C e~ 2lkT uT 2 (\ — u 2 ) 1+2lk , or equivalently where 
e~ 2jkT <C (1 — u) 1 ~ 2jk = (1 — u) 2 l k . So our approximation in the intermediate 
region is valid where 1 — u 3> e~ klkT . Expansion (3.7) in the parabolic region 
shows that 1 — u e~ kjkT = e~( fc /2-i)T precisely when (7>1. 

Going in the other direction, as u \ 0, one has the expansion 

(3.18) Q„[Zx] a -2c 2 (n - A) U - 6 + Ac 2 {n - 3)(1 + 2 7fe )^ 4 

- 2c 2 (l + 2 7fe )[n(l + 4 7fc ) - 10 7fc - l]u~ 2 . 
Therefore, in general dimensions, the solution Z 2 of equation (3.13) satisfies 

(3.19) Z 2 {u) = 0{u- i ) as u\0. 
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Near the tip, where u rs 0, the first two terms in our formal expansion of z predict 
that 

(3.20) z = c k e-^ kT u- 2 + 0{e- i ->" T u~ i ). 

This shows that our approximation in the intermediate region is expected to be 
valid in general only for u 3> e~ 7fcT , or equivalently for -0 3> e~( 1 ~ 1 / fe ) T , which is in 
agreement with equation (6.11) of [2]. 

3.3. Approximate solutions at the tip. Since the solutions in the intermediate 
region are predicted to be valid only for u ^> e~ lhT , we introduce the variable 



to be used in the tip region, which is defined in equation (5.14) below. Then, in 
view of (3.11), one has 



dz 

and hence 



dz 



- lkTZ r 



dz 
dr 

We rewrite this as 1 r [z] — 0, where 



2(n-l) £ " M " k Zr = 2(n-l) £rW " k Zr 



r 1 I 



— E r [z]. 



(3.21) 7 r [z] , . , „ , k .., } 2{n _ iy 

For large r, we assume for the formal argument here that we can ignore the terms 
containing e~ 2 ~ tkT 1 which reduces the equation to 

£ r [z] = 0. 

The solutions of this equation (see [1, 2]) are the Bryant solitons, 

z(r) = <&{a k r), 

with a k > an arbitrary constant (to be fixed below by matching considerations). 
To obtain more accurate approximate solutions of T r [z] = 0, one could try an 
expansion of the form 

(3.22) z = <8(a k r) + e" 2 ^ 7 ~P 1 (r) + e~^- T f3 2 (r) + ■■■ . 

In Section 5.2, we construct sub- and supersolutions based on the first two terms 
of this expansion. If as in [1] and [2] we normalize the Bryant soliton by requiring 

(3.23) Q5( r ) = l±|ill ) (r ->«>), 

then the Ansatz (3.22) leads to 

z ~ a A T 2 r -2 (r — > oo), 
and thus in terms of the u coordinate, to 

z ~ a A 7 2 e" 27fcT u~ 2 (r — > oo, u small). 
Matching this expression with that given in (3.20), we obtain c k = a k 2 , namely 

(3.24) a k :=^—. 



DEGENERATE NECKPINCHES IN RICCI FLOW 



11 



3.4. General features of the formal solutions. In addition to recalling (with 
certain adaptations) the explicit coordinate representations of the formal solutions 
in the parabolic, intermediate, and tip regions, we wish to note a few general 
features of the formal solutions that play a role in the statement and proof of our 
main theorem. We first note that the matched asymptotics analysis of [2] holds for 
all space dimensions n + 1 > 3, for all Hermite indices k > 3, and for all negative 
real values of the constant bk- We thus observe that the set of formal solutions is 
parameterized by these three numbers {n, k, bk}- 

Next, it follows from the form of the formal solutions in the parabolic region 
(and from the properties of the Hermite polynomials) that the formal solutions are 
reflection-symmetric across the (s = 0) equator if and only if k is even. These 
are the solutions that have degenerate neckpinches forming simultaneously at both 
poles. Also, asymptotically, they have only tip, intermediate, and parabolic regions. 
Hence, in studying (below) the convergence of full Ricci flow solutions to these 
formal approximations, it is only for the solutions with odd k that we must work 
in an outer region (beyond the parabolic region) as well as in the other regions. 

Finally, we note that if we calculate the curvature for one of these formal solutions 
(see Section 6 of [2]), we find that the norm of the curvature tensor achieves its 
maximum value at the tip, where it takes the value 

(3.25) | Rm(tip)t) | = _^L_ 

with Ck a constant 3 depending only on n, k, and bk- The curvature of a Ricci flow 
solution that asymptotically approaches this formal solution necessarily blows up 
at the same rate. 

4. Outline of the proof 

The main result of this paper, as stated above in the main theorem, is that for 
each formal degenerate neckpinch solution described above (with specified values of 
the parameters {n, k, bk}), there exist rotationally symmetric Ricci flow solutions 
that asymptotically approach this solution, and therefore share its asymptotic prop- 
erties. This result verifies the conjecture in Section 7 of [2]. The rest of this paper 
is devoted to proving this result. In this section, we outline the overall strategy we 
use to carry out the proof. In the remaining sections, we provide the details. 

Our strategy for proving that there exist solutions obeying the asymptotic pro- 
files derived in [2] follows Wazewski's retraction principle [13], as used in [5]. We 
apply this to solutions of Ricci flow, regarded as trajectories 



g : [0,T ) -> C°° ((T*S' i+1 ® Bym T*S" +1 ) 



'+> 

satisfying the isometry condition T*g(t) = g(t) for all T £ SO(n+l) and t £ [0, To). 
As noted in Section 2, such solutions may be naturally identified with ordered pairs 
of maps 

(<p,il>) : [0,T )-> (R+ xR + ) 
satisfying the boundary conditions (2.4) for all t £ [0, Tq), where To = Tq(<?o) is the 
maximal existence time of the unique solution with initial data go. 



3 At the tip, the sectional curvatures 3C T of a 2-plane tangent to an S n factor and X 1 - of an 
orthogonal 2-plane both take the value X = — — ■ with Ck '■= fc2 ( n 4 _ 1 )2 (2|fefc |)~ ¥ ■ 
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To implement this strategy, we define tubular neighborhoods S e of the formal 
(approximate) solutions discussed above and in [2]. For a given fixed formal solution 
9{n,k,b k }(t) with singularity at time T, and for a given positive e <C 1, we choose 
a tubular neighborhood of g{ n ,k,b k }(t) by specifying time-dependent inequalities to 
be satisfied by tp and ip, with these inequalities becoming more restrictive as t /* T, 
narrowing to the formal solution at time t = T. Decomposing the boundary of the 
tubular neighborhood 92 e = (<9S £ )_ U (dE £ ) U (dE £ ) + , we construct barriers and 
prove entrapment lemmas which show that solutions g £ S e never contact (9E e )_|_. 
We also prove an exit lemma which shows that any solution g £ S e that contacts 
must immediately exit. This ensures that the exit times of those solutions 
that never reach the neutral part (9S £ ) of the boundary depend continuously on 
the initial data. In Section 7, we construct (for each specified formal solution) a k- 
dimensional family of initial data g a (to) £ (S e PI {t — to}) parameterized by a £ 23 fe 
(the closure of a fc-dimensional topological ball) such that all g a (to) with a £ <923 k 
lie in the exit set (9H e )_, and such that the map : a £ d"B k —> g a (to) € (<9S 6 )_ 
is not contractible. (As seen in Section 7, fo £ [0, T) is chosen to be very close to 
the singularity time T of the formal solution.) We choose the initial data g a (to) 
at a sufficiently large distance from the neutral boundary (9S £ ) to guarantee that 
as long as the corresponding solution {g a (t) : t > t } stays within E e , it never 
reaches (9S e ) . In particular, the g a {ta) are chosen so that the solutions g a (t) 
can only leave S £ through (<9S e )_. If for every a £ 23 fc the solution g Q (i) were 
to leave S e at some time t a > to, then the exit map a i-> g a (t a ) would yield a 
continuous map $ : 23 fe — > (<9S e )_. Since t a — to for all a £ d% k , the map $ 
would extend and therefore ^ would have to be contractible, in contradiction 
with our construction of <\>. It follows that for at least one a £ 23 , the solution 
g a (t) never leaves the tubular neighborhood S e . Such a solution must then exhibit 
the asymptotic behavior derived in [2] (and summarized above in Section 3) for the 
formal model solution g{ n ,k,b k }(t)- 

We outline further the steps of this strategy in Section 4.4 below. Before doing 
so, we briefly describe how the tubular neighborhood S e (hereafter called "the 
tube") is defined in each of the four regions. In subsequent sections, we show that 
the tube, as defined here and (more explicitly) in Section 5, has all the properties 
needed to carry through the proof. 

4.1. Defining the tube in the tip and intermediate regions. In both of 
these regions, we define the tubular neighborhood S E using upper and lower barrier 
functions, z + (u,t) and Z-[u, r), which are defined with respect to the coordinate 
system described in Section 2.2. By the parabolic maximum principle, an evolving 
solution cannot cross a barrier. Thus the barriers provide an entrapment condition, 
ensuring that solutions never exit through the portions of (<9S e ) + they define. The 
barriers are constructed for a precise subset of u £ [0,1) using functions adapted to 
each region separately. To complete the construction, we "patch" these functions 
together in the transitional areas where the regions overlap, using the fact that the 
maximum (minimum) of two subsolutions (supersolutions) is a subsolution (super- 
solution). This work is done in Section 5. 

4.2. Denning the tube in the parabolic region. If S e could be defined by 
barriers everywhere, it would follow that solutions starting from an open set of 
rotationally symmetric initial data would all develop degenerate neckpinches at the 
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Figure 1 . A schematic view of a tubular neighborhood S e . 



it « 1 + b k e- x " T h k (a) 




Figure 2. An overview of the separate regions in the solution 



same time t — T. This is of course impossible, since one can always change the 
singularity time of a solution by shifting it in time. So rather than constructing 
barriers in the parabolic region, we must adopt an alternate (necessarily more 
complicated) approach. 

We define S e in the parabolic region working with the coordinate system dis- 
cussed in Section 2.1. In order to do so, we must introduce more notation. As 
discussed in Section 3.1 and in [2], the operator A defined in (3.4) naturally arises 
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in the analysis of the parabolic region, and we are led to write the formal solution as 
an expansion in its Hermite polynomial eigenfunctions. To make this more precise, 

2 

we note that A is self-adjoint in the weighted Hilbert space \) := L 2 (R; e _T da). 
Denoting the inner product in f) by (/, g)^ := J R f(a)g(a) da, and the norm 
by ||/||[) ■— %/(/, /)[,, we find that the Hermite polynomials {h m (-)}^ =0 constitute 
a complete orthogonal basis of f). We normalize so that h m (a) — a m + • •• . 

Solutions u = u-ldo not belong to fj, because they are not defined for all a. So 
we fix r\ > to be a smooth even cutoff function such that r)(y) — 1 for < \y\ < 1 
and rj(y) = for \y\ > |. We then define 



(4.1) v(a,r) 



r n(e-f kT / 4 a) v(a, r) for \a\ < §e 7fcT / 5 , 
forlcrl > |e 7feT/5 . 



It follows that v € f) is smooth and satisfies v(a,r) — v{a,r) for \a\ < e lkT ^ 5 . One 
computes from equation (3.2) that 



(4.2) 
where 



= Av + i]N[v] + E[ij, v), 



lrr , , s 2(n - l)v 2 - v 2 rr n 
N[v](a,r) := n/ 

#[77, u](ct,t) := (?/ r - 77^ + ^77^ v - 2r) a v 



with 



Vt = -^ae-^ r / 4 77'(e-^ T/4 (7) 



and where I[v] = I[u] is defined in (3.3). We note that in calculating r]N[v] and 
E[rj,v\, one may take v to be identically zero for any \a\ > |e 7fcT / 5 at which v is 
not defined. 

Because v e f), it makes sense to define the projections 
(4.3a) - fc -:=g W ^ 
(4.3b) Ufc:= ^_^fr fc , 

(4.3c) v k+ := 2^ fej- (= u - (u fc _ + Ufc)) . 

With regard to the linear operator A, the projections Ufe, and v k + represent 
rapidly-growing, neutral, and rapidly-decaying perturbations, respectively, of the 
formal solution. We show below that Vk- is associated to the boundary (9S e )_, V). 
to the boundary (dS E ) , and Vk+ to the boundary (<9S e ) + . 
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In the parabolic region, s e is denned by four conditions: a solution g(-) belongs 
to S e if it satisfies the three L 2 inequalities 

(4.4a) \\v k -h<ee-^ T , 

(4.4b) \\v k -b k e~^ T h k \\ h <ee-^ T , 

(4.4c) \\v k+ U<ee-^ T , 

together with the pointwise inequality 

(4.4d) sup (M + \d a v\) < We~^ T , 

\ a \<p 

where b k is the parameter that selects a particular formal solution, P is a constant 
chosen below to define the range of the parabolic region, W is a large constant to be 
fixed below, and [i k is defined in (3.5). We show below that the first inequality above 
is an exit condition, meaning that any solution that contacts <9S e with = 
£e -fj. k T immediately exits 3 e . We further show that the remaining three inequalities 
are entrapment conditions, ensuring that solutions of interest never contact the 
parts of 9S £ that they define. This work is done in Section 6. 

4.3. Defining the tube in the outer region. If k > 4 is even, we do not have 
to deal with an outer region when defining S £ : in this case, the formal solution 
described in Section 3 is reflection symmetric. So we encase it (outside the central 
parabolic region) in barriers that are themselves reflection symmetric, placing such 
properly ordered and patched barriers on either side of the parabolic region. 4 

If k > 3 is odd, on the other hand, we define S e outside the parabolic region with 
respect to the coordinate system described in Section 2.2, using barriers adapted 
to a precise subset of u € (l,e T / 2 ). This work, which is highly analogous to the 
constructions for the tip and intermediate regions, is also done in Section 5. 

4.4. Summary of the main estimates used in the proof. We presume now 
that we have fixed n > 2 and k > 3, and chosen b k < 0, thereby selecting a formal 
solution u(cr, r) = 1 + b k e^ llkT h k (a) + ■ • • (corresponding to g{ n , k .b k } ( r )) to serve as 
a model in the parabolic region for the solutions whose existence we prove in this 
paper. Two parameters remain to be determined: the small constant e appearing 
in (4.4a)-(4.4c), and the large constant W appearing in (4.4d). Carefully tracking 
which estimates depend on which choices, we show below that at the end of the 
proof, we can choose W sufficiently large, depending only on b kl such that our exit 
and entrapment results, Lemmas 15-18, and hence the Main Theorem, hold for all 
sufficiently small e, depending on {b k , W}. Toward this end, we proceed as follows. 

Using the relation z = i\) 2 = 2(n— l)u 2 , we employ properly ordered and patched 
barriers constructed in Section 5 to argue in Section 6 that there exist exist con- 
stants A , B, 5i, and Tg, 5 depending only on b k , such that solutions v = u — 1 in 
the tube S e satisfy 

^<|a CT {( e «>|)i}|<^ 



^If k > 4 is even, it is not necessary that the initial data we construct in Section 7 be reflection 
symmetric. For such data to belong to H s , it suffices that they lie inside the reflection symmetric 
barriers defined in Section 5, and that they satisfy estimates (4.4a)-(4.4d). 

^At several steps of the proofs in Sections 5—6, one has to restrict to sufficiently large times t; 
we address this in Section 7 by constructing initial data that become singular sufficiently quickly. 
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for t > t s and \a\ £ ( 



|P,(5ie 7fcT ), where 



(4.5) 




) 



k 



The inequality above provides a precise sense in which solutions in S e resemble 
the formal solution for large \a\, i.e., in the interface between the parabolic and 
intermediate regions. Accordingly, we define the parabolic region as 



For later use, we note that P and the other constants introduced in Section 5 
depend only on bk and are in particular independent of W. 

Further analyzing the parabolic-intermediate interface, we prove in Section 6 
that there exist C± depending on {bk, W} such that solutions in E £ satisfy 



in y for t > Tg, where C± depends only on bk, and r g > t 8 depends on {bk,W}. 
We next argue that there exists C2 depending only on bk such that one has an L 2 
bound 



for all r > t 10 , if T10 > t 9 is chosen sufficiently large, depending on {bk, W}. 

To finish the proof, we apply our estimates for ||(9 T — to the three L? 

inequalities (4.4a)-(4.4c) defining S e in the parabolic region to prove an exit lemma 
for the first inequality and entrapment lemmas for the remaining two. Then we use 
the analytic semigroup generated by A to derive C 1 bounds 



that hold at all sufficiently large times r > t\\ > T10, where t\\ depends on {bk, W}, 
but C3 depends only on bk- We use this estimate to prove an entrapment lemma 
that says that for W large enough and all small enough e, solutions in S e never 
contact dE e by achieving equality in (4.4d). Collectively, these results hold for 
all times r > f > th, thereby establishing all necessary properties of the tubular 
neighborhood S e . The details of the steps outlined here for the parabolic region 
are provided in Section 6. 

The argument we have outlined above proves that for at least one initial datum 
belonging to S e at the initial time, the solution that flows from this datum develops 
a degenerate neckpinch — provided that the set of such initial data is nonempty. 
The work to construct such data is done in Section 7, where we choose to £ [0, T) 
close enough to T so that t = to corresponds to r = f. This completes the outline 
of the proof of the main theorem, which validates the predictions we made in [2]. 6 
We proceed in the sections below to supply the details of the proof. 



One may translate between the coordinate representations used in [2] and those used here by 
following the methods illustrated in Section 3 above. 



(4.6) 



y := {a : \a\ < P}. 




\\(d T -A)v\U<C 2 e-^ 



sup {\i(a,T)\ + \i a (a,T)\} <C 3 e-» 
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5. Constructing barriers 

5.1. Barriers in the intermediate region. In the intermediate region, one says 
that z is a subsolution (supersolution) of equation (3.11) if D u z < (> 0), where 



!/ -1 \ -2 Q u[z] 

— {u — u)z u — u z • 



(5-1) V v z:= Z , .. 

or u 2 2(n - 1) 

A lower (upper) barrier is a subsolution (supersolution) that lies below (above) a 
formal solution in an appropriate space-time region. These are obtained by first 
constructing parameterized families of sub- and supersolutions, and then making 
suitable parameter choices to ensure properly ordered barriers. 

Lemma 1. Let Z\{u) = Aiu~ 2 (l — u 2 ) 1+2lk for some A\ > 0. There exist a 
function £ : (0, 1) — ¥ K, a constant B± > 0, and a constant A™ ln < oo depending 
only on A\ such that for any A2 > A 1 ™ 11 , the function 

(5.2) z ± (u, T ) := e' 2 ^ Z x {u) ± A 2 e~ i ~"' T '£(«) 

is a supersolution (+) or a subsolution (— ) in the space-time region 

fa~ / a \ k i 2 

<M ____ „ , n / A 2 



(5.3) Bi y^ e ~ lkT - u - 1 - Bi ( ^ ) e ^ fcT ' T - Tl ' 

where T\ depends only on A\ and A 2 . 

We do not claim that C( u ) > 0, so the subsolutions and supersolutions provided 
by this lemma may not be ordered. In Lemma 2, below, we construct properly 
ordered barriers z± by specifying a pair of constants A^ and A^ , setting A\ = A^ 
in the expression (5.2) for z+ and A\ — A^ in the expression (5.2) for z_. The 
time-dependent domain (5.12) where these barriers are valid provides the precise 
characterization of the intermediate region, as promised in Section 3.2 above. 

Proof. We verify here that z + is a supersolution. (The verification that z_ is a sub- 
solution follows from the same argument, with the sign of A 2 changed.) Applying 
the operator D u defined in (5.1) to the expression (5.2), we obtain 



[z+] = e-^j-^u" 1 - u)Zj! - {u- 2 + 2 lk )Z l ) 



+ e- 47fc7 



A2 -e-^Q u[ZuC] _ _H_ e -47^ Qti[c] |. 
1 2(n — 1) J 



2(n 

It follows from the definition of Z\ that the terms multiplying e~ 2lkT vanish; thus 
we have 

(5.4) e 4 ^V u [z + ] = A^- 1 -^- 1 -u)C - (u' 2 +4 7fe )c} - ^Ti) ^ 

An A 2 

_ A2 e -27 fc rp ft ^2 



n-l 2(n-l) 

We want to choose C and A 2 so that for sufficiently large values of r, CD u [z + ] > 
(2)„[z_] < for the subsolution). To achieve this, we first note that as r — > 00, 
only the first two terms on the right in (5.4) remain. We therefore choose £ so that 
those first two terms together are positive so long as A 2 is large enough. This then 
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implies that for any fixed u £ (0, 1), one has D u [z+](u, r) > if r is large enough. 
We conclude the proof by estimating how large r must be for the constant terms 
in (5.4) to outweigh the remaining terms. 

Estimate for Q u [Zi]: For < u < 1 (a necessary condition for the intermediate 
region), it is easy to see that 

C C 

\Zi'\ < —r kcZi and \Z X "\ < ^Z x . 

u(l — u- i ) u /i (l — u z ) z 

Using Q u [Zi] = ZiZ x " - \{Zi) 2 ~ vr 1 Z l Zi - 2(n - l)u~ 2 Zf, we find that 

(5.5) \Qu[Zi]\ < CA\u-\l - u 2 ) 4 ^. 

Here and in what follows, unnumbered constants C are generic and may change 
from line to line. 

Choice of (: Let £ : (0, 1) — > K be any solution of the inhomogeneous ode 

(5.6) - \{u^ - «)^ - (4 7fe + u- 2 )C = u- 6 (l - u 2 )^. 
Replacing u by u 2 as the independent variable, one can rewrite this equation as 

v ' A(u 2 ) I 1 - u 2 u 2 J u 

whereupon it is easy to see that the equation has regular singular points both at 
u = and at u = 1. The solution of the corresponding homogeneous equation is 
given by (3.14) with m — 2, namely Chom = c 2 uT 2 (\ — u 2 ) 1+4lk . One can either 
use the Frobenius method of power series or variation of constants to conclude that 
every choice of C solving (5.6) has the asymptotic behavior 

u- 4 + O(u- 2 \ogu) (u\0), 
k -(1 - u 2 ) 4 ~ lk + 0((1 - u 2 )^ k+1 log(l - u 2 )) (u /■ 1) 

at the endpoints of the interval < u < 1. 



(5.8) CM 



Constructing super solutions: Estimate (5.5) and formula (5.4) together tell us that 
e 4 ~* kTr D u [z + ] > {A 2 - CA\}u- & {\ - u) 4 ~ tk 

-e-^^Z^QW - ^^e-^\Q u [C}\. 



n-1 1, ^ u 2(n-l) 

To ensure that the r-independent term is positive, we must choose A 2 > CA\. We 
set 

Af n = 2CA 2 , 
which then implies that for all A 2 > A™ 111 , one has 

(5.9) e 4 ^ T V u [z + }> 

CA 2 lU -\l~u) 4 ^ - J^. e -a7*r|^^ c]) _ -^i-e-^HQnlCll. 
n — 1 2(n — 1) 
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Determining where z + is a supersolution: By definition, one has 

2Q U [Z U C] = ZiC + Z/'C - Z/C - u-^Z.C + Zi'C) - 4(n - l)u" 2 ZiC. 
For < u < 5, we recall (5.8) to see that 

\Qu[Zi,C]\ <CA lU ~ 8 

and 

|Qu[C]| <Cu- w . 

Combining these estimates, we find there is c > such that 

e^ kTr D u [z+] > cu- & (A\ - CA 1 A 2 e- 2 ' 1kT u- 2 - C A 2 e^< kT u^) 



cu 



-G 



.4 



2 fl - C M - C ^ ^ 
1 V AAei^uY A 2 (e^uY> 



A 1 (e^ T u) 2 A\(e^ kT u) 4 

for < u < g. It follows from this inequality that there exists a constant B\ < 00 
such that z + is a supersolution in the region (e 7fcT u) 2 > B 2 A 2 /A U equivalently 

Bl ^_ e <.<_. 

On the other hand, for the case | < tt < 1, we have Z\ ~ Ai(l — u 2 ) 1+27fc and 
C - (1 - m 2 ) 4 ^, so that 

|Q„[^i,C]| <C^i(l-u 2 )- 1+6 ^ 

and 

iq«[c]| <c(i-7i 2 r 2+8 ^. 

Thus in the interval g < w < 1, we have 

c(l - u 2 )^ k {A\ - CA 1 A 2 e' 2 ^ kT (l - u 2 )- 1+2 ^ k - CA 2 e~^ kT (1 - w 2 )- 2+4 ^}. 
Hence, making B\ larger if necessary, we see that z+ is a supersolution if 

1/(1-27*) / a \ fc /2 

Since 1 — u 2 = (1 — it)(l + u), we have 1 — u < 1 — u 2 < 2(1 — w), which proves the 
lemma for z + . The same arguments apply to z_, so the proof is complete. □ 

We now construct a pair of properly ordered barriers, based on the sub- and 
supersolutions of Lemma 1. To do this, we first specify 

(5.10) Af:=(l±6)c k 

where is defined in (3.10), and < S -C \ is a small constant fixed here once 
and for all. We then define ^3 := max{A2(Aj~), A 2 (Af)}, where the dependence 
of A 2 on A\ and A\ is that which is discussed in Lemma 1, with the added 
requirement that ^2(^4^) > A 2 nm (Af ) are large enough to satisfy condition (5.23) 
in Lemma 5 below. (This is necessary so that the ordered barriers we construct 
here are compatible with those we construct in Section 5.2 for the tip region.) Wc 
further fix r 2 := max{T 1 (A^),T 1 (Af )} and B 2 := max{B 1 (A~), B 1 (A+)}. We then 
have the following: 



20 



SIGURD B. ANGENENT, JAMES ISENBERG, AND DAN KNOPF 



Lemma 2. There exist r 3 > ti and B3 > B2 depending only on b k , such that if £ 
is the function appearing in Lemma 1, then 

(5.11) z ± (u,t) = (1 ± (5)c fc e- 27fcT M- 2 (l - w 2 ) 1+27fc ± A 3 e-^ kT C{u) 
are a pair of properly ordered barriers in the space-time region 

i<2 a (1A \ k l 2 

(5.12) B2 \j~ - e ~ lkT < u < 1 ~ B 3 ( — -J e"^, T>r 3 . 

Proof. Using the asymptotic expansion (5.8) for £, we easily determine that z_ < z + 
as u \ 0. It is also straightforward to verify that on any interval c < u < 1 — c, 
there exists r 3 > T2 such that z_ < z + if c < u < 1 — c and r > r 3 . 

To study the rest of the intermediate region, we define w := 1 — u 2 . Then as 
u 1, the asymptotic expansion (5.8) implies that 

w" 1 " 27fc (z+ - z_) > 25c k u- 2 - 2A 3 e- 4 ~ tkT w 2 ~ lk ~ 1 (l + wlogw). 

Using the facts that 2j k — 1 = — 2/fc and that \w log w| < 1 for < u> < 1, one then 
estimates for u < 1 — _B 3 (2A 3 /cfe) fe / 2 e^^' i ' r and r > r 3 that 

w" 1 - 27fc (z + - z_) > 26c k u~ 2 - AA 3 e~^ kT w- 2/k 

> 2c k {5u- 2 - B- 2/k (l + u)- 2 / k e^ kT3 } 

> 

provided that B$ > B2 is chosen sufficiently large. □ 

5.2. Barriers in the tip region. In Section 3.3, we observe that the evolution 
of the metric is governed by the parabolic PDE 7 r [z] = 0, where r — e 7fcT it, and 
where 7 r is defined in equation (3.21). This motivates us to look for sub- and 
supersolutions — and then properly ordered upper and lower barriers — of the 
form z± — < B(A4r)±e^ 2 ' fkT f3(r), where /3(r) is to be chosen to solve a suitable ODE, 
where A4 is a parameter to be determined, and where we recall that 58 denotes the 
functional expression for the Bryant soliton. 

Lemma 3. For any A4 > ; there exist a bounded function j3 : (0, 00) — > R, a 
sufficiently small B± > ; and a sufficiently large T4 < 00, all depending only on 
A4, such that 

(5.13) z± := <B(A 4 r) ± e~ 27fcr /3(r) 
are sub- (z~) and super- (z+) solutions in the region 

(5.14) < r < B A e lkT 
for all t > T4 . 

We note that inequality (5.14) provides a working definition of the tip region, 
which overlaps the intermediate region (5.12) for all sufficiently large times t. 

As with the intermediate region, we do not claim that j3 has a sign; to get a pair 
of properly ordered barriers for the tip region, one makes suitable choices of the 
parameter A 4 (namely A^ and A\ as defined in (5.20) below), relying on the fact 
that 03 is monotone decreasing. 
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Proof. Using the abbreviation B(r) := 93(Air), we see that for the function z = 
B(r) + e~ 2lkT (3(r) to be a a supersolution, it suffices that 7 r [z] > 0. Using the 
definition of 7 r in (3.21) and the definitions of L r and Q r from (2.8) and (2.9), 
respectively, we find that 

(5.15) T r [B(r) + e~ 27feT /3(r)] = 

^ 2lkT {\rB'{r) —±—(L r[ p{r)] + 2Q r [B(r), P(r)})} 



2(n 

< T {-2 lk P{r) + r -f3'(r) - ^-L_Q r [/3(r),/?(r)]}. 



Here we use the fact that B satisfies 8, r [B] = 0. For large r, the sign of the RHS 
is determined by the sign of the coefficient of e~ 2 " <kT . This coefficient is linear 
inhomogeneous in j3 and its derivatives. To ensure that B(r) + e - 2 ^ kT /3(r) is a 
supersolution for large enough r, we choose f3(r) in such a way that the factor 
multiplying e~ 27fcT is positive. More precisely, we recall that rB'{r) < for all 
r > 0, which leads us to choose (3 to be a solution of the linear inhomogeneous 
differential equation 

(5.16) L r [(3(r)} + 2Q r [B(r),p(r)} = 2(n - l)ArB'(r), 

with the constant A to be chosen below. Applying the definitions of L r and Q r , we 
find that (5.16) takes the following form: 

(5.17) %) d » + {^-EV)-^}« 

&r z I r r ) dr 

r * R'(r) 1 — 2R(r) ~i - - 

+ [B"(r) y~ + 2(n - 1) 7 ^ Ai }/?W = 2(n - l)ArS'(r). 

This equation has a two-parameter family of solutions. We seek a specific solution 
from this family that leads to the positivity of the RHS of (5.15) for sufficiently 
large r. To determine the region on which B(r) + e~ 2lhT (3(r) is a supersolution, we 
need to know the asymptotic behavior of our solution /3(r) both as r \, and as 
r Z' co. We now show that this asymptotic behavior is the same for all solutions 
of (5.17), because near both r = and r — oo, we find that the solutions to the 
homogeneous equation are much smaller than any particular solution. 

Asymptotic behavior at r = 0: At r = 0, we have B(0) = 1, B'(0) — and 
B"(0) = -^«8"(0). Since 93 (r) is an analytic function at r = 0, the ODE (5.17) 
has a regular singular point at r = 0, and to leading order can be written as 

d 2 l(r) + n^2 d#r) _ 2(n - 1) = ^ 
dr 2 r dr r 2 

for some constant C. The general solution of this equation is 

Mr) = cur 1 - 11 + a 2 r 2 - Cr\ 

for some C. Discarding the unbounded solution and choosing a 2 — 1, we conclude 
that there exists a solution fi p (r) of the true ODE (5.17) that satisfies p (r) — 
r 2 + o(r 2 ) as r \ 0. 
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Asymptotic behavior atr = oo; We have normalized 03 so that Q3(r) = r~ 2 + 0(r~ 4 ) 
as r — > oo. This choice yields B(r) = ^S(A^r) = A± 2 r~ 2 + 0(A^ 4 r~ 4 ) near r = oo. 
Keeping only the lowest-order terms in r , we see that for r large, the ODE (5.17) 
is a small perturbation of the ODE 

1 I. . , n - 1 „, , , 2(n - 1) „ . , 4(n - 1)A 
(5.18) ]4p/&(»0 + — + -^T-t-Pooir) = K , 



(5.19) ^ 0o (r) = a 1 re- ar "+a 2 r / p^e^ >dp- ^ 2 . 



whose general solution is 

12' 

M 

where a := ^^A 2 . The first two terms in (5.19), which solve the homogeneous 
equation associated to (5.18), both vanish as r / oo: the first has Gaussian decay; 
the second is 0(r~ 2 ). Hence every solution of the true ODE (5.17), including the 
choice of (3 p (r) made above, has the property that j3(r) = (1 + o(l))( — 2A/A 2 } as 
r —> oo. 

To determine how to choose A, we now estimate 



%[B{r) + e~ 2lkT p p {r)] = e~ 2 ~ (kT (\ - A)rB'(r) 



+ e- 4 ^ T (-2 7fe /3 p (r) + ^/3;(r 



Q r \P P (r),P p (r) 
2(n- 1) 



Since B'(r) < 0, the first term on the RHS is positive if we choose A > jr. To make 
a definite choice, we set 

Near r = 0, oo, the asymptotics of 93 then imply that 

_ £, n = jV 2 + o(r 2 ), (r\0) 
r l7 " J \c 2 r 2 + o(l/r 2 ) (r^oo) 

for certain positive constants ci and c 2 , and hence that 

—rB'(r) > cminjr 2 , ^ 2 }, 

for some constant c. We also know for j3 p that 

2A 

/3 p (r) = 0(r 2 ) (r\0), /? p ( r ) ^ - + (r" 2 ) (r / oo). 
These expansions can be differentiated. Hence we have for < r < 1, 



-2%/3p(r) + r#(r) 



r fl/ ^ Q r [^(r),/3„(r)] 



< CY 2 



; 2(n- 1) 
and thus 

7 r [B(r) + e- 2 ^ kT l3 p (r)} > - e - 2lkT rB \r) - e~^ kT Cr 2 
> e- 27fcT r 2 (5-e- 27feT C) 
>0, 
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if t is large enough. For r > 1, we have 

Q_r/=LM 

<C, 



7fePplj+ fc^ lj 2(n-l) 



so that 

T r [B(r) + e- 27feT /3 p (r)] > -e" 27tr rB'(r) - e" 47fcT C 

> e -27 fe r( 6r -2 _ C r e -2 7 »Tj 

>o, 

provided that r < <5e 7fcT , where J = \Jc/C. 

This completes the construction of the supersolution. A subsolution can be 
constructed along the same lines, but one could also just set /3(r) = and observe 
that (5.15) then takes the form 

7 r [B(r) + e- 2 ^ T p{r)] = \e~ 2 ^ r r& {r). 

Because B'(r) < for all r, the RHS of this equation is negative; hence we see that 
B(r) — 03 (^Ur) is a subsolution. □ 

We now use the sub- and supersolutions from Lemma 3 to construct properly 
ordered upper and lower barriers for the tip region. To do this, we specify the 



constants Af as follows, using the small constant < S -C 1 fixed above: 



(5.20a) (A4 )- 2 = {1 + S)(l+ ^B^ (1 - 6)c k , 

(5.20b) (A+)- 2 = (1 - S) (l + l -B^ (1 + 8)c k . 

We note several consequences of these choices. First, they ensure that both barriers 
are close to the formal solution Q3(c fc r). Second, they ensure that A~\ < A±; 
because 03 is monotone decreasing, this implies that *B(A^r) < *B(A%r) for all 
r > 0, which is useful in establishing the desired ordering of the barriers. Third, 
the choices above, and in particular the leading 1±5 factors, ensure that A A satisfy 
the strict inequalities in condition (5.22) of Lemma 5 below. We use this fact to 
guarantee that the barriers constructed here are compatible with those constructed 
in Lemma 2 within the intersection of the intermediate and tip regions. These 
regions overlap for all sufficiently large r, depending on Bi in Lemma 2 and B4 in 
Lemma 3. 

Lemma 4. If the parameters A 4 take the values specified by (5.20), and if (3 is the 
function appearing in Lemma 3, then there exists T5 > T4 such that 

(5.21) z ± = 03(A±r) ± e~ 2lkT (3{r) 

form a pair of properly ordered barriers that are valid for all points < r < B±e lkT 
and all times r > r 5 . 

Proof. Near r = 0, one has *B(A^r) — 1 — b(A^ r) 2 + 0(r 4 ) for some constant b > 0, 
and j3(r) = [1 + o(l)]r 2 . Thus as r \ 0, we see that 

z+ - z_ = {b[(A- A ) 2 - (A+) 2 } + 2[1 + o(l)]e- 7 ^}r 2 + 0(r 4 ) > 0. 
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Near r = oo, one has the expansions *B(Afr) = 1 + (Af r)~ 2 + 0(r~ 4 ) and 
[3 = [1 + o(l)]{ - 2(1 + l/k)(A±y 2 }. It then follows from the specifications of Af 
given in (5.20) that 



1 - 2[1 + o(l)]^e- 2 ^ r | + 0(r- 4 ), 



which is positive for all sufficiently large r and r. 

It is straightforward to verify that z + > z_ on any bounded interval c < r < C, 
for all sufficiently large r. The result follows. □ 

5.3. The intermediate-tip interface. Let 2™* denote the upper and lower bar- 
riers constructed in Lemma 2 for the intermediate region, and let z± p denote the 
upper and lower barriers constructed in Lemma 4 for the tip region. The domains 
of definition of these two pairs of functions intersect for sufficiently large times. To 
obtain barriers that work throughout the union of these regions, we must verify that 
zlr?* and z± p are properly ordered in their intersection. In what follows, we focus 
on upper barriers, omitting the entirely analogous argument for lower barriers. To 
avoid notational prolixity, we write Ai and A4 for Af and A± , respectively, in the 
statement of the lemma and its proof. 

Lemma 5. Suppose that A\ and A4 satisfy strict inequalities 

(5.22) (1 + Ib^A, < A-, 2 < (1 + \b^)A 1} 

where B 3 is the constant from Lemma 2. 
Then there exists a constant C such that if 

(5.23) A z >C^A\ /2 + {l+ l -Bf) 2 A\ 

then one has A 3 > C(A\^ 2 + A^ 4 ), and consequently 

(5.24a) 4 P < atr = B 3 ^A 3 /A 1 , 

(5.24b) z+ p > z^> atr = 2B 3y /A 3 /A 1 , 

for all t > tq, where Tg > max{r3, T5} is sufficiently large. 

Proof. Near infinity, the Bryant soliton has an expansion of the form (see Lemma 18 
of [1]) 

B( P ) = i f + 4+*§ + ... 
p 2 P A p 6 

Written in terms of the r coordinate, the supersolutions satisfy the expansion 

zl p = <B(A 4 r) + e^ kT (3{r) 

= Al 2 r- 2 + b 2 A± A r- 4 + 0(r- 6 ) + 0(g- 2 ^ T ) 

as r — > 00. The 0(r -6 ) term comes from the asymptotic expansion of < S(A 2 r) and 
is therefore uniform in time (in fact, independent of time). 
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Using our asymptotic description (5.8) of £, we find that as u \ 0, one has 

Z™* = e- 2 "< kT Z x {u) + e- A ~< kT A z C,{u) 

= A ie - 2 ~< kT u- 2 (l - u 2 ) 1+2 ~< k + A 3 e- 47fcT i*- 4 (l + 0(u 2 logu)) 

= Air~ 2 + A 3 r-* + 0(re- 27fcT ), 

where r — e 7fcT u, and where 0(re _27t:T ) is uniform on any compact r interval. 
Hence on bounded r-intervals, one has 

r 2 (zf - zf) = (Af - A,) + (b 2 A^ 4 - A 3 )r- 2 + 0( r - 4 ) + 0(re-^ kT ) 



= (Al 2 - A,) + (b 2 A± 4 + 0(r- 2 ) - A 3 )r- 2 + 0(n 



We want this quantity to change from negative to positive as r increases from 
B 3V /A 3 /A 1 to 2B 3y /A 3 jA 1 . At r = B a ^A a /Ai, we have 



a: 



Ax 



b 2 A! 

A 3 A\ 







AAB. 



0(re" 



and at r = 2B 3 a/ A 3 ~[A\ , we have 



2 / 
r [z_ 



tip 
+ 



v int\ 



A: 



A, + 1 I b2Al 











^1x2 

aJ 



A 1 \ S 3 - 2 + 0(re~ 27fcT ). 



If A 3 > CA 4 and A 3 > Cy/Ax both hold for sufficiently large C, then 



A 3 Ai +U[ A 2 > 



< 



Ai 
2 ' 



and wc find that 



r 2 {zf 



2/ tip 



< A 



> A 



4 

-2 



(1 + ^ 2 )A, + 0(re- 2 ^) at r = B 3 ^A 3 7T 1 , 
(1 + '^B^A, + 0(re- 2 ^ T ) at r = 2B 3 ^A 3 ~[A[. 



From this it is clear that if the strict inequalities (5.22) hold, then the upper barriers 

□ 



zV p and zT 1 satisfy the patching condition (5.24) for all sufficiently large r. 

5.4. Barriers in the outer region. If k > 4 is even, then as explained in 
Section 4.3, solutions in 5 e are controlled by the barriers constructed above. If 
k > 3 is odd, we need additional barriers to define (<9S £ ) + in a suitable subset of 
ue (l,e T / 2 ). 

Lemma 6. For any A$ > 0, there exist a function £ : (l,oo) — > K and constants 
Aq and B$ > depending only on A§ such that 

(5.25) z± = A b e- 2 ^ kT {u 2 - l) 1+27fc U - 2 ± A 6 e~^ kT C{u) 
are sub- (z"—) and super- (z + ) solutions in the region 

(5.26) 1 + B 5 e-^ kT < u < B^e T/2 . 

Proof. We sketch the proof, because it is very similar to that of Lemma 1. 

Then for 1 < u < oo, one estimates \zi\ < 



Let 5i(u) = (u 2 - lY+^i 
C{u-l)^ k+1 u^ k '\ \z[\ < C{u-l) 2 ~* k u 2 ~< k -\ and \z'{\ < C{u - l) 2 ™" V^" 1 . 
Hence 

\Qu[h]\ <C{u-\)^ k u^ k ~ 2 . 
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This motivates us to choose the function £ to be a solution of the ODE 

h u _ u -i)^ _ (4 7fc + U - 2 )C = (u- If^u**- 2 . 
2 au 

Its solutions are £ = Cp + Zi. h om , where £ p is any particular solution, and i^hom 
solves the associated homogeneous ODE. We may thus choose a solution £ such that 
C ~ (u — l) 4 T fc (the behavior of Cp) as it \ 1 and £ ~ w 87fc (the behavior of ^2, hom) 
as it -> 00. Estimating that |(j < C(w - l) 4 ^u 4 ^, |C'| < C{u - l) 4 ^" 1 ^, and 
IC"I < C(u - l) 4 T fc - 2 u 4 ^, we observe that 

\Q u [Q\<C(u-l) 8 ^- 2 u 8 ^, 

and 

|Q«[2i,C]| < C( u - l) 6 ^- 1 ^" 1 . 

With these estimates for Zi(tt) established, we write 

2± = ^ 5 e- 27fcT zi(u) + A 6 e~ 4 ^C( U ), 

as above, and seek to determine Aq (depending only on A$) so that z + is a super- 
solution. Based on the estimates above, we have 

This leads us to choose A e := 1 + 2 (n-i) ^i' Then one obtains 

e 4 ^^„[z + ] > (u - l) 4 ^7i 4 ^- 2 {l - CA 5 A 6 X - C(A 6 X) 2 }, 

where X := er 2lkT (u — l) 2 ~ fk ~ 1 u 2 ~ <k+1 . It follows that z + is a subsolution in any 
spacetime region where CA 5 (A e X) + C(A & X) 2 < 1. 

For u > 1 near 1, we observe that there is a constant B' depending only A§ such 
that z + is a supersolution if (w - 1) 1_2 T* > B' e - 27fcT , hence if w 2 /' £ > B' e - 27fcr . It 
is easy to see that this is implied by the constraint u > 1 + Bc,e~^ kT if B5 is large 
enough. 

For u > 1, we observe that there is a constant B" depending only on A$ such 
that iL|_ is a supersolution if B"u < e r / 2 , hence if < 1. 

Similar considerations show that z_ is a subsolution. □ 

As in the other regions, by making suitable modifications to the constants, one 
can ensure that the sub- and supersolutions of Lemma 6 are properly ordered upper 
and lower barriers in the outer region. To do so, we specify constants Af = (l±<5)cfc, 
where is defined in (3.10), and < 5 <C 1 is the small constant fixed above. We 
define A 7 := max{A 6 (ylg ), A 6 (Af)}, and we set B 6 := max{B 5 (A^),B 5 (A^)}. 
Then we have the following: 

Lemma 7. There exist B 7 > B e depending only on such that if £ is the function 
appearing in Lemma 6, then 

(5.27) z± = (l± 5)c fe e- 27fcT (u 2 - l) 1+2 ~< k u- 2 ± A 7 e-^ kT ((u) 
are a pair of properly ordered barriers in the region 

(5.28) 1 + B 5 e-^ kT <u < (B 7 y 1 e T/2 . 
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Proof. We sketch the proof, which is similar to that of Lemma 2, but somewhat 
simpler. 

As it \ 1, one calculates that £ ~ [u — l) 47fc > 0. And as u — > oo, one has 

e -2 7fc r^2 _ 1 jl+2 7fctt -2 _ e -2 7fc r w 4 7fc and e ~4 7fc r^ _ {eT 2 ^ U^») 2 . It follows that 

the barriers are properly ordered provided that B-je~ T l 2 u < 1 for some B-j > Bq 
sufficiently large. □ 

6. Analysis of the parabolic region 

Ricci flow solutions cannot escape from the portion of 9S e that is associated with 
the tip, intermediate, and outer regions, and is explicitly constructed in Section 5 
using upper and lower barriers. Solutions can, however, escape from the portion of 
9S e that is associated with the parabolic region, and is defined by the inequalities 
listed in (4.4a)-(4.4d). In this section, after establishing a collection of key estimates 
in Lemmas 8-14, we derive the exit and entrapment results, Lemmas 15-18, that 
relate to these inequalities and control S £ in the parabolic region. As noted in 
Section 4.4, these results play a crucial role in our proof of our main theorem. 

6.1. The parabolic-intermediate interface. We recall from equation (3.9) that 
in the parabolic region, the formal solution indexed by n, k, and bk satisfies 

z « 2(n - l)fc 2 (-& fc ) 2/fc e~ 27fcT (l - u) 1+2 ^ 

for | | 3> 1; here we use the identity 2 — 2/k = 1 + 2jk. As shown in Lemmas 1-5, 
we can construct properly ordered and patched barriers that encase this formal 
solution and are valid for r > tq and u < 1 — B%e~^ kT \ i.e., for e Mfcr |'i;| > Bg. 
Here tq and B$ := B^{2A^/ cu) k ^ 2 depend only on bk- These barriers define S e in 
the intermediate and tip regions. We now use them to derive information about 
the parabolic- intermediate interface for solutions belonging to S e . This interface 
corresponds to |cr| sufficiently large, where "sufficiently large" is made precise below. 

Lemma 8. There exist constants A s , B% and t-j > tq depending only on bk such 
that every solution in S e satisfies 

\l A~~ \ A + 

(6.1) < d a {(-e^v)*} < y-j± 

for B 8 < e' IfeT H < \e^ kT and r > t-j. 

Proof. By using (5.8) to bound |£| < C(l — u) Alk for 4 < u < 1, we see from 
Lemmas 1-2 and the fact that 27^. = /i^(l — 27^) that there exist constants A~, A + 
such that any solution in S e satisfies 

A-e- 2 ^ T u- 2 (l - u 2 ) 1+2 -< k < z < A+e- 2 ^ T u- 2 (l - u 2 ) l+2 ^ 

for i < u < 1 — B$e ~^ fcT and r > t-j > t$, where A , B$, and t-j depend only on 
bk ■ This implies that there exist constants Ag depending only on bk such that one 
can bound = fc- 2 (l — u) * 2 z in the same space-time region as follows: 

fe~^ <(^)) 2 <|K 2 ^ 

Recalling that — v = 1 — u > and (ik/k = 7&, we obtain (6.1). □ 
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We now show that Lemma 8 gives pointwise control of e Pfcr |w| for large |er|, 
specihcally in the parabolic-intermediate intersection. Below, we use properties 
of the operator A to provide pointwise control of this quantity for smaller \a\. 
Both of these arguments are facilitated by a judicious choice of the parameter 
P = 2(B /\bk\Y^ k defined in equation (4.5). Noting that there is a universal constant 
Cfc with \h k (a) — c | < C k \a\ 2 for |er| > 1, we choose B > B% large enough, 
depending only on b kl so that 

(6.2) P >^^c~ k . 

Lemma 9. IfP = 2(B/\b k \) 1 / k is chosen as in (6.2), depending only on b k , then 
for all sufficiently small e, depending on {b k ,W}, every solution in S £ satisfies 

d lkT \v(^P,T)\ < B, and e^ kT \v(a,r)\ > B if \p<\cf\<P. 

Proof. We first show that the bounds are true for the formal solution. Our choice 
of P ensures that 




On the other hand, if \a\ > |P, then 

\b k \h k (a) > \b k \\a\ k {1 - C k \a\- 2 } > \\b k \\a\ k . 

In particular, 

w4r)>(^B>l B . 

Now because II - II CT 2 < II • Ik, one may apply a Gagliardo-Nirenberg 

L 2 (9;e 4 dcr) 

interpolation inequality (e.g., see p. 125 of [11]) in the compact parabolic region 7 
to see that 

\\e^ T v-b k h k \\ L ^ {y) 

< c[\\d a {e^ T v - b k h k )\\)!l m \\e^ T v - b k h k \\ 2 / 3 + \\e^ T v - b k h k %] . 

Conditions (4.4a)-(4.4c) imply that \\e^ kT v — b k h k ^h < 3e for solutions in H e . Hence 
by condition (4.4d), one has 

||e"* T v - b k h k \\ t ~ m < C(l + W 1 / 3 )e 2 / 3 . 

So the result holds for all sufficiently small e, depending on {b k ,W}, where the 
constant W is chosen in Lemma 18 below. □ 

Remark. Because v is continuous, it follows that there exists tTs(r) g (f-P) f-P) 
such that e IJ ' kT v(aB(T), t) = —B, 

Hereafter, we assume that e — s(b k ,W) is sufficiently small, as indicated in 
Lemma 9. One then has the following result, anticipated in Section 4.4: 

Lemma 10. There exist = A 8 , 5±, and t& > t-j depending only on b k such that 
for t > rs, every solution in S e satisfies 
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in the parabolic-intermediate interface \a\ € (|P, 5ie 7fcr ), where P = 2{B /\bk\) 1 ^ k ■ 

Proof. Lemma 9 implies that any solution «£H E obeys the derivative estimate of 
Lemma 8 for all |er| > |P such that \v\ < \. Further, so long as Lemma 8 applies, 
one has 

(e»i> T \v\)i < B.i + (\a\ - ~P* 



k V 5 

i 

The RHS is bounded above by (|e /ifcT ) k provided that \a\ < 8ie lkT and t > r 8 , 
where we choose r§ > T7 large enough so that 8\e lkTa > 2P. □ 

We next establish the following bound, also anticipated in Section 4.4: 

Lemma 11. There exist C± depending on {b^, W} such that for \o~\ £ (§P, 5ie 7fcT ) 
and t > t$, every solution in !B e satisfies 

C_\a\ k < e"* r |v| < C+\n\ k . 

Proof. With C_ := ( min{S 1 /' £ /(|P), J /k}) k , we apply Lemmas 9-10 to see 

that (e^ kT \v\) 1/k > C Vfe da for §P < a < 5ie^ kTl . An analogous estimate holds 
for -<5ie 7fcri < a < -f P in the case that k is even. This proves the lower bound. 
For the upper bound, note that (e Alfer |w(|P, t) |) 1 / fe < B x l k holds by Lemma 9. 

In the region |P < |er| < |P, one has \bk\hk(o-) > ^ (|) fe P for the formal solution, 
and hence e^ T \v\ > (§) + B, as in the proof of Lemma 9. So for |P < |er| < |P, 
condition (4.4d) implies that d a (e^ kT \v\) h < where 

If we now choose C+ := ( max{B 1 / fc /(§ P), T?(6 fc ), /fc}) , then the upper bound 
follows. □ 

Remark. Combining Lemma 10 with Lemma 11, one obtains 

(6.4) C'Jo-f- 1 < e^ T \v a \ < C' + \a\ k - 1 

for \a\ G (|P,(5ie^ r ), where C' ± := Cl~ 1/k VA± depend on {b k ,W}. 

6.2. Estimating the "error terms". The results in the previous section prepare 
us to derive useful bounds for \(d T — A)v\ and \\{d T — 

Lemma 12. There exist C\ depending only on b k , and tq > t% depending on 
W} 7 such that for every solution in S E , the pointwise bound 

\{d T -A)v\ < dW 2 e- 2flkT 

holds in the parabolic region T for all t > tq . 

Proof. Recalling that v = v in the parabolic region, we determine from (4.2) that 
(d T — A)v = N\ oc — nv a I in T, where 

2(n- l)v 2 a - v 2 



loc 



2(l+i;) 
and 

7= / T \ da. 

1 + v(cr, t) 
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If we choose r 9 > r 8 , depending on {bk,W}, so that We~^ kT9 < \, then it 
follows from (4.4d) that one has 1 + v > 1 — We~~^ kT > \ for \a\ < P and r > r 9 . 
Inequality (4.4d) also implies that there exists Co (depending only on bk) such that 
|Moc| < C W 2 e~ 2f * kT holds in 7 for r > r 9 . 

Integration by parts shows that 



(6.5) / = 



1 



,2 



1 



da. 



In 7, the first term on the RHS of (6.5) is bounded by 4We ^ kT , while the second 
is bounded by 2PW 2 e~ 2fJ - kT < APWe~^ kT . Noting that \v a \ < We~^ kT in 7, we 
obtain a suitable bound for \v a I\ and thus complete the proof. □ 

Lemma 13. There exist Ci depending only on bk, and t\q > tq depending on 
{bk, W}, such that for all t > t\q, every solution in S e satisfies 



\\{d T -A)v%<C 2 e 

Proof. It follows from Lemma 12 that there exists Co depending only on bk such that 
a pointwise bound \{d T — A)v\ < Coe~2 A " tT holds for \a\ < P, and for all sufficiently 
large r, depending on {bk, W}. Thus to complete the proof, we concentrate on the 
rest of the support of (d T — A)v, namely P < \o~\ < |e 7fcT / 5 . We assume that r is 
large enough so that e 7fcT / 5 < <5 1 e 7fcT . 

It follows from equation (4.2) that (d T — A)v — r/(N loc — nv a I) + E for \a\ e 
(P, |e 7fcT / 5 ), where E = (rj T — r\ aa + f^) v - 2r) a v a . 

For e chosen as in Lemma 9, it follows from Lemma 10 that 1 + v > | in the 
interval (|P, i5ie 7fcT ) . In the same region, one has C_|cr| fc < e Atfcr |'i;| < C+|fj| fc and 
C'_\cr\ k ~ 1 < e^ kT \v a \ < C' + \o-\ k ~ 1 as consequences of Lemma 11 and estimate (6.4), 
respectively. 

Combining these inequalities, and again integrating by parts to evaluate /, we 
obtain C\ depending on {bk, W} and C2 depending only on bk such that one has 

\N loc \ + \v a I\ < C ie -^ kT \a\ 2k < C 2 e-^ kT 

in the interval P < \a\ < |e 7fcT / 5 , at all r sufficiently large, depending on {bk, W}. 
Here we use the fact that 2^7^/5 < /itfc/2. 

Similarly, the pointwise estimates for v and v a above further imply that there 
exist constants C3,Ci depending on {bk,W}, and C5 depending only on bk, such 
that the estimates 

/>oo 

\\E\\l < C 3 e- 2 ^ T / \o-\ 2k+2 e-° 2 /*da 

/ „27 fc r/5\ 

< C 4 e- 2 ^ T exp f - B -^ r J e (»+i)7*r/5 

< Cse" 3 ''^ 

hold for all t sufficiently large, depending on {bk, W}. The result follows. □ 

Finally, we derive pointwise bounds for v that are independent of W — bounds 
that apply to solutions originating from initial data we construct in Section 7 below. 
Given a smooth function f(o~) and a constant R > 0, we define 

WfWcHR) ■■= sup (\f(a)\ + \U(a)\). 

\tr\<R 
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Lemma 14. If at time t%i > txq, one has 

\\H;r 11 )\\cH2P) <Me~^ 
for some constant M > 0, then for all r > tu, one has 

(6.6) \H;t)\\ c i {p) <C 6 (l + M)e~^ T , 

where Cq > 1 depends only on (and is in particular independent ofW). 

Proof. By (4.2), we have v T — Av = F(a,r), where F — r]N[v] + E[r],vj. Thus for 
t > Tu + 1, the variation of constants formula lets us write 

(6.7) v(;r)=e A v{-,T-l)+ f T e (T - T ' ) - A F(-,r')dr'. 

Jt-1 

Standard regularizing estimates for the operator A imply that 

||e^(-,T-l)|| cl(P) <C 1 \\v(;r-l)\\ t) , 
and for < r — r' < 1, that 

|| e (r-r ')*F{W)\\ci(P)<&(T-T'r 3/ *\\F(W)h, 

where C\ depends only on P, hence by (4.5), only on bk- Conditions (4.4a)-(4.4c) 
imply that there exists C 2 depending only on bk such that 

||«(-,r-l)|| 6 <C 2 e~^ T . 

Furthermore, because rn > rio = Tio(bk,W), Lemma 13 provides C2 depending 
only on bk such that 

p7(-,T / )|| 6 <C a e-i'* r/ . 

Combining these estimates and integrating (6.7) gives (6.6) for times r > t\\ + 1. 
To prove (6.6) for t € (t 11; tu + 1), we again use variation of constants, writing 

(6.8) fi(-, r) = e (r - Tll) ' /t w(-, r u ) + / e {r - T ' )A F{- , r') At'. 

J Tu 

As above, we can use the estimate for \\F(-, t')||j given by Lemma 13 to obtain a 
satisfactory C 1 estimate for the integral term. However, we cannot use the smooth- 
ing properties of the operator e eA to estimate the first term on the RHS, since the 
time delay 9 = t—t\\ may now be arbitrarily short. But since A is a nondegenerate 
parabolic operator, there exists C3 depending only on P, hence only on 6^, such 
that for 9 € (0, 1), one has 

l|e M w(-,Til)|| C l (P) < tf 3 ||iK-,Til)|| C l(2F). 

Since \\v(-, Tn)||c ,1 (2P) < Afe _AlfcTl1 by hypothesis, this lets us estimate the first 
term on the RHS of (6.8). Thus we obtain (6.6) for m < r < ru + 1, which 
completes the proof. □ 

6.3. Exit and entrapment results. We now use the estimates derived above 
to prove the exit and entrapment results corresponding to various portions of the 
boundary of 5 e in the parabolic region. Our first result shows that solutions im- 
mediately exit if they contact (t?H s )_. 

Lemma 15. There exists f > t 10 depending on {bk,W} such that if a solution 
t)£S E contacts <9S e by achieving equality in (4.4a) at r > f , then it immediately 
exits H e . 
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Proof. We recall from Section 4 that the projection Vk- represents the rapidly- 
growing perturbations of a formal solution. Since the projection from v to 
commutes with the operator A, and since the projections {i>fe_, Ufe, Vk+} are pairwise 
orthogonal, one has 

' " / ' h ' ) = {v k -,Av+{d T -A)v\. 

a) h 



We readily verify the inequality {vk-^Av)^ > — fik-i\\vk-\\f,, from which it follows 
that 

iA|| e ^ fe _||2 > (Mfc _ ^iie^r^ua _ | (^v/» T (ar-A)«), |. 

Hence using Lemma 13 and Cauchy-Schwarz, we obtain 
1 d 



2^11^-118 



>e(--C»e-^ T ) >0 

e"* T t;,,_|U=e V 2 ' 



at all times t > r, for t > rio chosen sufficiently large, depending on {bk, W}. □ 

Our next three results are entrapment lemmas. 

Lemma 16. There exists f > t\q depending on {bk, W} such that solutions v £ S e 
cannot contact 9S e 6y achieving equality in condition (4.4c) ai any r > f. 

Proof. Arguing as in the proof of Lemma 15, one obtains 

Here, one has (vk+jAv)^ < — from which it follows that 

l__d 
2dr' 

Combining this inequality with Lemma 13, we see that if a solution were to contact 
9S e , then one would have 



:\K kT v k+ \\l < (p k - ^ + i)\\e^ T v k+ \\l + I (e^v k+ ,e^ T (d T - A)v\ |. 



I — Wp^t ||2 

2dr l|e vu+U 



<e(^-+C 2 e-^ kT ) <0 



for all r > t > t\q. By continuity, in a neighborhood of {||e AIfcr ?;/ c -|_||f, = e}, one has 
^r||e Alfcr Vfc+||^ < for t > f; this ensures that equality is never achieved. □ 

Our final two entrapment results apply to solutions in originating in the set _DJ 
defined by 



(6.9) Dl :-- 



IM-,T)-& fe e-^ fe ||o < |e^ fc " 
||ffe(-,T)-6 fc e-^^ fc || c i(2P) < 100e"« f 



This set is designed so that solutions satisfying v(-, r) £ S e D D T e cannot exit S £ by 
violating conditions (4.4b) or (4.4d). 

Lemma 17. There exists f > t 10 large enough, depending on {bk,W}, such that 
solutions satisfying v(-,f) £ S e PI D r e cannot contact <9S e by achieving equality in 
condition (4.4b) at any t > f. 
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Proof. We first note that one can write v k (a, r) = 4>{t)t){P l a)h k (cr), where the 
quantity 

<Kt) :=IIMI^ 2 (fi,fc fc ) e 

can be shown to evolve by 

^<t> = -tl k <t> + \\h k \\ h 2 ((Or - A)v, h k \ . 

Observing that 

\\e^ T v k ~b k h k U = (e" T #-&fe)||/i fc ||(,, 

one computes that 

l^W^'vk - b k h k \\l = \\h k \\l(e^ T 4> ~ h^i^t-h) 

< vWe^Vk - hhkh\ {(d T -A)v,h k \ \. 
It then follows from Lemma 13 that one has 

^-\\e^ T v k - b k h k \\ h < C 2 \\h k \\t,e-^" T . 
clr 

So if v £ S £ n D T e , then for all r > f , one has 

\\e^v k -b k h k \\,^+ 2 ^ke- 3 ^ <e 

provided that f is chosen sufficiently large, depending on b k and e (which depends on 
{b k , W} via Lemma 9). We conclude that equality in (4.4b) cannot be attained. □ 

Lemma 18. There exists W sufficiently large, depending only on b k , and there 
exists f > Tio sufficiently large, depending on {b k , W}, such that solutions satisfying 
«(-, t) £ S £ fl DJ cannot contact <9S e by achieving equality in (4.4d) at any r > f . 

Proof. Applying Lemma 14 with m = f > t\q and M = 100 yields 

sup {\v(a,r)\ + \v„(<t,t)\} < C 3 e-^ T 
W\<p 

for all r > f, where C3 := 101C6 depends only on b k . Hence the result holds for 
any W > C 3 . □ 

7. Constructing suitable initial data 

With the properties of the tubular neighborhood S e surrounding a given for- 
mal solution g{n tk ,b k }(t) established (primarily in Sections 5-6), including the exit 
and entrapment properties of portions of the boundary of H e , the remaining work 
needed in order to complete the proof of our main theorem is to show that there 
exist parameterized sets of initial data yielding Ricci flow solutions that have the 
properties discussed in the introduction to Section 4. In particular, presuming that 
the choices of 9{ n ,k,b k }{t) have been fixed, we seek a continuous bijective map $ 
from a closed topological fc-ball 23 fe to the set of rotationally symmetric metrics 
on S' i+1 with the following properties, for some choice of time to £ [0, T): (i) the 
image of $ is contained in S e p|{t = to}; (h) the image of $ restricted to d% k is 
contained in the exit set (cG £ )_; (iii) ^la^fe is n °t contractible; and (iv) the image 
of $ is far enough away from the neutral part of the boundary (9S £ ) to guarantee 
that so long as a Ricci flow solution that starts from initial data in $(23 fe ) stays in 
S e , it never reaches (<9S £ ) D . 
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We show in this section that for each choice of gi n ,k,b k } (*) anc ^ corresponding 3 e , 
sets of initial data satisfying these conditions can be found, if e is small enough. 

We first fix the initial time to at which we choose our sets of initial data. Noting 
that the results stated in the exit and entrapment lemmas discussed in Section 6.3 
are only guaranteed to work for r > f (corresponding to t > T — e~ r ), we set 
to := T — e _T , where f may be presumed to be large enough so that one has to > 
and | e T^/ 5 > IP. 

5 — 

We next define initial data locally using coordinates adapted to the parabolic 
region of the formal solution, since it is in this region where the restrictions stated 
above are most critical. Working with the metric variables and coordinates dis- 
cussed in Section 2.1, and recalling that k > 3 and bk < are fixed, we see that for 
every set of k constants 

bje[-B r B 3 ], 0<j<k-l, 

(with Bj and Bj small positive constants to be determined below), if we write 

k - 

u = 1 + E bje-^h k (a), (V| < ~e^/ 5 J , 

3=0 

and if in addition we set f) = 1, then a metric g = {da) 2 + u 2 g can is determined 
locally. It is clear that one can choose 13 j ,Bj > small enough such that the 
corresponding metrics are as close as desired in the C 1 topology to the approximate 
solution l + bke~ t * kT hk(cr) in the region \a\ < |e 7fcT / 5 . Closeness to g{ n ,k,b k } in that 
same region then follows. 

With an eye toward conditions (4.4a)-(4.4c) and (6.9), we choose ]3j,Bj > so 
that the sharp inequalities 

H-he-^hkWt, < 

||w fc +|| 6 <ee-^ f , 
hold for (6o, .. . ,6*_i) G n^o ( - M.j,Bj), with 

||« fc _||„ =ee~^ 

holding on the boundary, which ensures by Lemma 15 that solutions originating 
from initial data on the boundary of the topological ball 11^=0 [ ~ —j'^j] imme- 
diately exit 3 e . We verify that in choosing B ■ and Bj as indicated, the conditions 
for these initial data to lie within the tube (at least for the portion of the metrics 
pertaining to the parabolic region) are satisfied. We may assume that e is small 
enough, hence that Bj, Bj > are small enough, so that the C 1 condition in (6.9) is 
satisfied. Combined with the second inequality above, this ensures that the initial 
data constructed here belong to D\. In turn, by Lemma 18, this ensures that the 
data satisfy (4.4d). 

We next obtain globally defined initial data in (S E [}{t — f }) |~| D T Z by using 
cutoff functions to smoothly glue the metric g = (der) 2 + u 2 g can defined locally 
above to the formal solution in the intermediate and tip regions, and (for k odd) 
to that portion of the manifold which remains nonsingular. More specifically, we 
proceed as follows. 
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For |e 7fcT / 5 < a < |e 7fcT / 5 , the construction above ensures that z = tp 2 is as 
close as desired to its value in the formal solution (3.9), namely 

HnMM^) = Cfce- 2 ^ T (1 - u) 1+2 ^ (1 > u » e" 7 ^). 

In particular, making e smaller if necessary, we can ensure that each initial datum 
lies between the barriers z_ < z + constructed in Lemma 2 for all u where both 
it and the barriers are defined. One may thus extend the solution by means of 
a cutoff function, obtaining z with the properties that z = z at M(|e 7feT/ ' 5 cr) and 
z = Z{n,fe,6 fc } I0r all < u < u(|e 7fcT / 5 cr). Here we use the fact that Z{ ni k,b k } 
smoothly extends to u = using the construction in Section 3.3. Moreover, Zi ni k,b k \ 
lies between the barriers constructed in Lemma 4. 

If k is even, we repeat this step for |e 7fcT / 5 < —a < |e 7fcT / 5 , and we are done. 

If k is odd and |e 7feT / 5 < — a < |e 7fcT / 5 , we again find that each initial datum 
z is as close as desired to the formal solution 

Z{n,k,b h }{u) = c k e- 2 ^ T (u 2 - l)u- 2 (l< u « e T ' 2 ), 

hence lies between the barriers z_ < z + constructed in Lemma 7. So we again 
extend z by means of a cutoff function until u — <5e r / 2 , namely tp = 5. There we 
smoothly glue on a smooth punctured sphere as in Section 8.2 of [3]. We omit 
further details. 
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